Production of very neutron-rich nuclei with a "^^Ge beam 
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Production cross sections for neutron-rich nuclei from the fragmentation of a ^^Ge beam at 
132 MeV/u were measured. The longitudinal momentum distributions of 34 neutron-rich isotopes 
of elements 13 < Z < 27 were scanned using a novel experimental approach of varying the target 
thickness. Production cross sections with beryllium and tungsten targets were determined for a large 
number of nuclei including 15 isotopes first observed in this work. These are the most neutron-rich 
nuchdes of the elements 17 < Z < 25 (^"CI, ^^'Ar, 55.56j^^ ^"^'^^Ca, 59.60,6ig^^ 62,63rpj^ 65,66y^ 68^^.^ 
^"Mn). A one-body Qg systematics is used to describe the production cross sections based on ther- 
mal evaporation from excited prefragments. Some of the fragments near ^*Ca show anomalously 
large production cross sections. 

PACS numbers: 25.70.Mn, 27.40.-l-z, 27.50.+e 
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I. INTRODUCTION 

The discovery of new nuclei in the proximity of the 
neutron dripline provides a benchmark for nuclear mass 
models, and hence for the understanding of the nuclear 
force and the creation of elements. Once the production 
methods and cross sections for the formation of neutron- 
rich nuclei are understood, further investigations to study 
the nuclei themselves, such as decay spectroscopy, may 
be planned. Therefore, the determination of production 
rates for the most exotic nuclei continues to be an im- 
portant part of the experimental program at existing and 
future rare-isotope facilities. 

A number of production mechanisms have been used 
to produce neutron-rich isotopes for 17 < Z < 25. For 
example, ^'^'^^V were observed for the first time in deep- 
inelastic reactions in 1980 Q, ^^^'^^V [I and ^^^V [| by 
fragmentation of a ^^Kr beam for the first time in 1985 
and 1992, respectively, and 62,63,64y ^qj-q g^g^ observed 
as products of projectile fission in 1997 4]. The on- 
line isotope separation technique was used to produce 
for the first time ^3^54 ^^^^ 53^^^ ^gg^ ^^^^ 
years later, 54,55,56Qg^ isotopes were observed following 
in-flight fission of ^ssy ^ -pjjg possibility of producing 
very neutron-rich isotopes (e.g. ^^'^°Ca) in transfer-type 
reactions at intermediate beam energies has been dis- 
cussed recently Q . 

New progress in the production of neutron-rich iso- 
topes is possible given the increased primary beam inten- 
sities at the National Superconducting Cyclotron Labora- 
tory (NSCL) at Michigan State University and advances 
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FIG. 1: (Color online) The region of the nuclear chart inves- 
tigated in the present work. The solid line shows the limit of 
bound nuclei from the KTUY mass model The 15 new 

isotopes observed for the first time in the present work are 
marked by red squares. 



in experimental techniques fi\ allowed probing further 
into unknown regions of the table of the isotopes with 
fragmentation reactions. Indeed, recent measurements at 
the NSCL 0, H, [§] have demonstrated that the fragmen- 
tation of ^^Ca and ^^Ge beams can be used to produce 
new isotopes in the proximity of the neutron dripline. 
Continuing this work, we report here the next step to- 
wards the fundamental goal of defining the absolute mass 
limit for chemical elements in the region of calcium. 

The search for new neutron-rich isotopes was carried 
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out with a primary beam of Ge and using the two-stage 
fragment separator technique that enabled the discovery 
of '^''Mg and '^^Al in 2007 In the present measure- 
ment, 15 neutron-rich isotopes with 33 < < 45 were 
identified for the first time (see FiglT]) and indications for 
the existence of a new island of inversion centered on ^^Ti 
— as predicted by Brown [ll[ — were described in a pre- 
vious brief presentation of this work [l^. In this paper, 
we describe the details of our experimental approach and 
discuss the results. 



for an additional discrimination against various scattered 
particles. 

The primary ^^Ge'^''+ beam current was continuously 
monitored by a BaF2 crystal mounted on a photomul- 
tiplier tube near the target position, as well as a non- 
intercepting capacitive pick-up before the target. These 
indirect monitors were normalized to Faraday cup read- 
ings throughout the course of the experiment. The av- 
erage beam intensity for the measurements of the most 
exotic fragments was approximately 1 e/xA. 



II. EXPERIMENTAL DETAILS 
A. Setup 

A 132 MeV/u ''^Ge beam, accelerated by the coupled 
cyclotrons at the NSCL, was used to irradiate a series of 
beryllium targets and a tungsten target, each placed at 
the target position of the A1900 fragment separator p^ . 
Previous studies Q demonstrated that using the A1900 
fragment separator as the sole separation stage signifi- 
cantly limits the sensitivity of the experimental setup for 
very rare isotopes due to 

• excessive rates and large pile-up probability in the 
detector located at the intermediate dispersive fo- 
cal plane of the A1900 fragment separator (required 
to measure the momentum), 

• high rates of light nuclei in the final focal plane 
detectors that cannot be rejected by the energy- 
loss technique. 

The combination of the A1900 fragment separator with 
the S800 analysis beam line Q| to form a two-stage sepa- 
rator system as described in Ref. 0] , allows a high degree 
of rejection of unwanted reaction products so that the full 
primary beam intensity can be utilized, while preserv- 
ing the unambiguous identification of single exotic ions. 
The first stage of the system serves as a selector whereas 
the second stage contains detectors and functions as an- 
alyzer, see Fig. O 

In the present work the momentum acceptance of the 
A1900 was varied from Ap/p = 0.1% to 5%. For all 
settings, the horizontal and vertical angular acceptances 
were ±65 and ±40 mrad. At the end of the S800 analy- 
sis beam line, the particles of interest were stopped in a 
telescope of eight silicon PIN diodes (50x50 mm^) with a 
total thickness of 8.0 mm. This provided multiple energy- 
loss measurements and thus a redundant determination 
of the nuclear charge of each fragment in addition to the 
total kinetic energy A 5 cm thick plastic scintillator po- 
sitioned behind the Si-telescope served as a veto detector 
against reactions in the Si-telescope and provided a mea- 
surement of the residual energy of lighter ions that were 
not stopped in the Si-telescope. A position sensitive par- 
allel plate avalanche counter (PPAC) was located in front 
of the Si-telescope. This position measurement allowed 



TABLE I: Flight paths between various timing detectors. 



N 


Timing signals between 


Length (m) 


1 


FP scintillator 


2"'' PIN detector 


46.0 


2 


OBJ scintillator 


3'''* PIN detector 


21.0 


3 


FP scintillator 


OBJ scintillator 


25.1 


4 


Target* 


3'''* PIN detector 


81.5 



* - from the arrival time relative 

to the phase of the cyclotron rf-signal 



The time of flight (TOF) of each particle that reached 
the detector stack was measured in four different ways 
(see Table|T|. All TOF signals were used for the rejection 
of unwanted events. The first three TOF measurements 
listed in the table were used to deduce the mass-to-charge 
ratio A/q with a better resolution than with a single TOF 
measurement. 

The magnetic rigidity {Bp) of each particle was ob- 
tained by combining position measurements from two 
PPACs located in the dispersive plane of the S800 anal- 
ysis beam line with NMR measurements of the dipole 
fields. 

The simultaneous measurement of multiple AE sig- 
nals, the magnetic rigidity, the total energy, and the 
TOFs of each particle provided an unambiguous iden- 
tification of the atomic number, charge state and mass 
number of each ion. The detection system and particle 
identification was calibrated with the primary beam and 
confirmed by the locations of gaps corresponding to un- 
bound nuclei in the particle identification spectrum. In 
addition, a germanium 7-ray detector was placed near the 
Si-telescope (see Figl2]) to provide an independent verifi- 
cation of the isotope identification via isomer tagging as 
described in Ref. llSl . 



B. Experimental runs 

The present experiment consisted of three parts that 
are summarized in Table |TT1 During all runs, the mag- 
netic rigidity of the last two dipoles was kept constant 
at 4.3 Tm while the production target thickness was 
changed to map the fragment momentum distributions. 
This approach greatly simplified the particle identifica- 
tion during the scans of the parallel momentum distribu- 
tions. 
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FIG. 2: (Color online) Sketch of the experimen- 
tal setup at the NSCL. This overview diagram 
shows the equipment used for the production, 
separation and detection of new neutron-rich 
nuclei. 



S800 analysis beam line 



TABLE IL Experimental settings 



Data 


Fragment 




Magnetic rigidity, Bp{Tm 


) 


Target 


Stripper 


Wedge 




Time 


Beam 


Goal 


set 


of interest 


D1D2 


D3D4 






DioDii 


mg/cm? 


mg/ cm} 


mg 1 cm?" 


(%) 


hour 


particles 




1 


43s 


4.3233 


4.3233 


4.3134 


4.3036 


4.3000 


Be 9.8 






0.1 


0.98 


1.41el4 




2 














Be 97.5 






0.1 


2.03 


1.50el4 


momentum 


3 














Be 191 






0.1 


1.98 


2.07el4 


distribution 


4 














Be 288 






0.1 


1.85 


1.17el4 




5 














Be 404 






0.1 


0.94 


5.76el3 




6 


43 s 


4.3372 


4.3233 


4.3134 


4.3036 


4.3000 


Be 191 




20.2 


1 


3.36 


1.25el5 


43m s 


7 














Be 288 




20.2 


2 


5.97 


1.96el5 


production 


8 


66y 


4.3609 


4.3383 


4.3222 


4.3059 


4.3000 


Be 404 




20.2 


5 


19.86 


1.42el6 




9 














W 567 


C 12 


20.2 


5 


20.74 


I.470I6 


production 


10 


54 Ar 


4.3529 


4.3332 


4.3192 


4.3051 


4.3000 


W 567 


Be 145 


20.2 


5 


10.33 


7.55el5 


of new 


11 














Be 629 




20.2 


5 


9.37 


6.15el5 


isotopes 


12 


59 Ca 


4.3566 


4.3355 


4.3205 


4.3054 


4.3000 


Be 629 




20.2 


5 


21.60 


1.58el6 




13 


58 Ca 


4.3546 


4.3342 


4.3198 


4.3052 


4.3000 


W 567 


Be 97.5 


20.2 


5 


46.59 


3.34el6 





The momentum acceptance of the A1900 fragment sep- 
arator was restricted to Ap/p = 0.1% for the first part of 
the experiment devoted to the measurement of differen- 
tial momentum distributions. The use of different beryl- 
lium target thicknesses (9.8, 97.5, 191, 288, 404 mg/cm^) 
allowed to cover different spans of the fragment momen- 
tum distributions necessary to extract production cross 
sections and also resulted in more isotopes in the par- 
ticle identification spectrum. The particle selection de- 
pends critically on the magnetic rigidity of the fragment 
separator because the average energy (and therefore the 
magnetic rigidity) of a given reaction product decreases 
with increasing target thickness. As a result, the isotopic 
selection of this system moved toward more neutron-rich 



fragments with increasing target thickness. 

For the second part of the experiment, a Kapton wedge 
with a thickness of 20.2 mg/cm^ (indicated as achromatic 
degrader in FiglJ) was used at the dispersive image of 
the A1900 to reject less exotic fragments with an 8 mm 
aperture in the focal plane while the separator was set 
for 43 s. The goal of this setting was to confirm the par- 
ticle identification by isomer tagging with 43m g ^j^^ _ 
319 keV, Ti/2 = 0.48 fis) and 32™A1 {E^ = 735 keV, 
Ti/2 = 0.2 fis). 

In the final part of the experiment, dedicated to the 
search for the new isotopes, four settings were used to 
cover the most neutron-rich isotopes with 13 < Z < 27, 
as it is impossible to find a single target thickness and 
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FIG. 3: (Color online) Particle identification plot showing the 
measured atomic number, Z, versus the calculated function 
A — 3q for the nuclei observed in this work {q is the ionic 
charge of the fragment). See text for details. The limit of 
previously observed nuclei is shown by the solid line as well 
as the locations of several reference nuclei. Particle identifi- 
cation was confirmed by registering gamma-rays from short- 
lived isomeric states of *''S and ^^Al (marked by black circles). 



magnetic rigidity to produce all of the fragments of in- 
terest. Each setting was characterized by a fragment on 
which the separator was tuned. A search for the most ex- 
otic nuclei in each setting was carried out with Be and W 
targets. The four settings were centered on ^^Ar, ^^-^^Ca 
and ^^V respectively, and were based on LISE++ [l^ cal- 
culations using the parameterizations of the momentum 
distributions obtained in the first part of experiment (see 
Section IIV Ap . The momentum acceptance of the A1900 
was set to the maximum of Ap/p = 5.0% for these pro- 
duction runs. In contrast to previous studies with the 
A1900 as single-stage separator, e.g. Q, the present over- 
all experimental efficiency was dominated by the accep- 
tance of the two-stage separator and not by the limita- 
tions due to high counting rates of unwanted nuclei. 



III. ANALYSIS OF EXPERIMENTAL DATA 

The advantage of our approach — keeping the mag- 
netic rigidity (Bpof) constant while varying the target 
thickness — can be seen in Fig. [31 which shows the total 
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(a) Q; a(Q) = 0.073 





(b) Z; a(Z) = 0.098 




FIG. 4: Ionic (top) and elemental (bottom) spectra ob- 
tained with Eq. IA3I and IA5I for all particles stopped in the 
Si-telescope during the data set dedicated to the search for 
new isotopes (see Table HI]). 



distribution of fully-stripped reaction products observed 
in all runs of this work. The range of fragments is shown 
as the measured atomic number, Z, plotted versus the 
calculated quantity A — 3g. The identification of the in- 
dividual isotopes in Fig. [3] was confirmed via isomer tag- 
ging using the known isomeric decays in ^^Al and '*'^S, 
and from holes in the A — = 1 line corresponding 
to the unbound nuclei ^^O and ^^F located between the 
particle bound nuclei ^■^N and ^^^Ne. The details of the 
calculation of the particle identification are given in the 
Appendix. 

The ionic (q) and elemental (Z) spectra obtained for 
particles that stopped in the Si-telescope are shown in 
Fig. 21 The peaks were fitted with Gaussian distributions 
of constant width. The standard deviations were found 
to be aq = 0.073 for the ionic spectrum and az = 0.098 
for the elemental spectrum, respectively. This means 
that the probabilities of one event being misidentified 
as a neighboring charge state or element are equal to 
PqiO.5) = 3.7 10-12 and Pz(0.5) = 1.7 10^^ respec- 
tively. 

The mass spectra for the isotopic chains from chlorine 
to manganese measured during the production runs are 
shown in Fig. O Only nuclei that stopped in the Si detec- 
tor stack are included in this analysis. The observed frag- 
ments include 15 new isotopes that are the most neutron- 
rich nuclides yet observed of elements 17 < Z < 25 

^50(^j 53^j. 55,56j^ 57,58Qg^ 59,60,61g(, 62,63rpj 65,66y 

^*Cr, ™Mn). The new neutron-rich nuclei observed in 
this work are those events to the right of the solid line 
in Fig. [3] and to the right of the vertical dashed lines 
in Fig. |5l As noted before [13], the previously reported 
observation of ^^Cl fT^ may have actually been an ob- 
servation of the hydrogen-like charge state ^^01^^+. 
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FIG. 5: (Color online) 
Mass spectra of the el- 
ements 17 < Z < 
25. All particles that 
were stopped in the Si- 
telescope during the pro- 
duction runs were ana- 
lyzed. The limits of pre- 
viously observed nuclei 
are shown by the vertical 
dashed lines. 
Standard deviations pro- 
duced in multi-peak fits 
with Gaussian distribu- 
tions at constant width 
(red dashed curves) are 
shown in the plots for 
each element. 



Mass number 



IV. RESULTS AND DISCUSSION 

A. Parallel momentum distributions 

It is important to be able to predict the momentum 
distributions of residues when searching for new isotopes 
in order to optimize the fragment separator at the max- 
imum production rate. Also, the accurate prediction 
of the momentum distributions allows estimation of the 
transmission and efficient rejection of strong contami- 
nants. Several studies of the parallel momentum dis- 
tributions have been made, e.g. [ll[ii,|23,IIl|, but the 
predictions for the production of the most exotic nuclei 
are still very uncertain. The few semiempirical mod- 
els [H, [l^ used to describe the data assume Gaussian 
momentum distributions, characterized just by two pa- 
rameters (the mean value and the width) that may not be 
sufficient to model the momentum distributions. There- 
fore, the measurement of the fragment momentum distri- 
butions remains an important part of the search for new 
isotopes. 

The most straightforward and common way to mea- 
sure the momentum distributions of projectile fragments 
is to scan the magnetic rigidity of the fragment separator. 
A thin target is generally used for this measurement to 
avoid complications from differential energy-loss in the 
target (the systematic change in the kinetic energy lost 
by the beam nucleus and product nuclei in the target). 
In addition, the momentum acceptance of the fragment 



separator is often restricted to a small value for isotopes 
with large yields. The acceptance is then opened up for 
the smaller yields, producing fewer measurements with 
larger systematic uncertainties. The total cross section 
is determined by integrating these momentum distribu- 
tions using some estimate, or measurement, of the angu- 
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FIG. 6: (Color online) Differential cross sections of "^Cr (top) 
and ^"K (bottom figure) fragments obtained by varying the 
target thickness at one magnetic rigidity. Red solid lines rep- 
resent the fitted Gaussian functions. Horizontal errors corre- 
spond to the velocity difference caused by production at the 
beginning or the end of target. 
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Thickness (mg/cm ) 

FIG. 7: (Color online) The measured yield of ^''Cr (top pan- 
els) and ^"K (bottom panels) as function of target thickness. 
The lines show the LISE"*"^ calculations with results of fit- 
ting the velocity distributions shown in Fig. |6] (green solid 
line) and default model settings [l^ (blue dashed line) as de- 
scribed in the text. Note that the left panels have linear hor- 
izontal scales, while the right panels have logarithmic scales 
to emphasize the difference between calculated curves for the 
thick and thin targets. 



lar transmission, e.g. j22l |. 

In the present work a new approach to measure mo- 
mentum distributions and cross sections was used. In 
contrast to the "Bp-scanning" method using one thin 
target, a variety of targets with different thicknesses was 
used at constant magnetic rigidity. This method is par- 
ticularly well suited to survey neutron-rich nuclei since 
the less exotic nuclei are produced with the highest yields 
and their momentum distributions can be measured with 
the thin targets. The lighter fragments with high yields 
will experience the largest differential energy losses and 
will fall outside the constant momentum acceptance of 
the separator with the thicker targets. Conversely, the 
heaviest fragments with the lowest yields will only be 
produced in sufficient numbers in the thickest targets and 
will have lower differential energy losses. Extensive simu- 
lations of the A1900/S800 tandem separator showed that 
this method can provide a sensitive measurement of the 
mean value and width of the momentum distribution [23j . 
However, as in the Bp scanning method, the momentum 
distributions of some fragments will be incomplete. 

Five targets were used to measure the momentum dis- 
tributions (see Table HI)) . The yields of two example frag- 
ments, ^^Cr and ^"K shown in Fig. [6] as a function of 
reduced velocity were fitted with Gaussian functions and 
the total production cross sections were inferred. The 
momentum distributions for 34 isotopes were derived (in- 
dicated by the bullets in Fig. [J) and integrated to deduce 
the cross sections. It should be noted that at the energy 
of these experiments, the shape of the fragment momen- 
tum distribution is asymmetric due to a low-energy ex- 



ponential tail stemming from dissipative processes [2lj. 
The Gaussian function used for the present data does 
not take this tail into account, but the underestimation 
of the cross section is small. As one example, the ^^Ca 
cross section is underestimated by approximately 2%. 

A survey of all of the fitted results showed that frag- 
ments in the heavy mass region were produced with 
significantly higher velocities and slightly broader mo- 
mentum distributions than the model predictions. The 
model [m used for these calculations assumes that the 
energy necessary to remove each nucleon is Es — 8 MeV, 
a value derived for fragments close to stability. Further 
analysis showed that the separation energy parameter 
for nuclei observed in the present work in the region 
Ap/2 < Ap < Ap exhibits a linear decrease with the 
number of removed nucleons: 

= 8- 11.2Ayl/^P (I) 

where — Ap — Ap^ Ap is the projectile mass number, 
and is the fragment mass number. 

The measured yields of the two example isotopes, 
^^Cr and ^"K, are shown as a function of target thick- 
ness in Fig. [T] The fragment yields calculated with 
the LISE++ code [l6] using parameters from Morrissey's 
model of the reaction with EFAX cross sections 24] are 
shown by the dashed lines. The fragment yields from 
the same calculations with the experimental parameters 
from fitting the velocity distributions (see Fig. [6]) are 
shown by the solid lines in Fig. [71 Only a fair agree- 
ment is found with the model calculations with default 
settings. On the other hand, the agreement between ex- 
perimental data and model calculations with parameters 
produced by fitting the velocity distributions is excellent 
and gives confidence that the fragment separator was set 
to optimize the yield of the weakest, most neutron-rich 
fragments. 

B. Production cross section 

The inclusive production cross sections for the ob- 
served fragments were calculated by correcting the mea- 
sured yields for the finite longitudinal and angular ac- 
ceptances of the separator system. A total of thirty-four 
cross sections with beryllium were obtained from Gaus- 
sian functions fitted to the longitudinal momentum dis- 
tributions; these nuclei are indicated by stars in Fig. [H 
The cross sections for all of the remaining fragments with 
incomplete longitudinal momentum distributions were 
obtained with estimated transmission corrections. The 
angular and longitudinal transmissions were calculated 
for each isotope in each setting using a model of the mo- 
mentum distribution with smoothly varying parameters 
extracted from the measured parallel momentum distri- 
butions. For example, the values of the angular, mo- 
mentum and overall transmissions, as well as the wedge 
selection for the central fragment in runs with a wedge 
are given in Table [1111 The wedge selection represents the 
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FIG. 8: (Color online) In- 
clusive production cross 
sections for fragments from 
the reaction of '^''Ge with 
beryllium and tungsten 
targets shown as a function 
of mass number. The cross 
sections with the beryUium 
targets derived by momen- 
tum distribution integration 
are shown by stars, those 
normalized with LISE^"*" 
transmission calculations 
are indicated by solid dia- 
monds The cross sections 
obtained with the tungsten 
target were normalized 
with LISE-^-*- transmission 
calculations. The two green 
dashed lines separate nuclei 
that require neutron pickup 
in the production mecha- 
nism. The red lines show 
the predictions of the EFAX 
systematics for beryllium 
(see text). 



Mass number 



fraction of isotopes that passed through the slits at the 
focal plane when the wedge is inserted at the intermedi- 
ate dispersive image. The estimate of the uncertainties in 
the transmissions came from the values at the one-sigma 
limits of the reduced width of the longitudinal momen- 
tum distributions scaled by the mass loss p^ : 

(7||=f7o%/A^ (2) 

where AA was defined above, and ctq = 105 ± 15 MeV/c. 
Details of transmission calculations, the analysis of their 
uncertainties, as well as the general analysis of momen- 
tum distributions taken from various target thicknesses 
will be described in a subsequent paper [23l |. 

The cross sections obtained for all of the fragments 



TABLE III; Transmissions during runs with a wedge 



Data 


Fragment 


Transmission(%) 


set 


of 


Momentum Angular 


Wedge 


Total 




interest 






selection 




6 


43g 


11 ±3 


70+^ 


7^-1-10 


6±2 


7 


43g 


16 ±3 


69^^ 


74^1° 


8±2 


8 
9 
10 


66y 
66y 

s^Ar 


70 ±3 
69 ±3 

47 ±6 


99+^ 
99tl 
88+t 


83tl, 
83+?3 

79l?4 


58l?o 
33tl 


11 
12 
13 


54 Ar 
59 Ca 
58 Ca 


30 ±3 
32 ±3 
411? 


8611 
93+1 
93l5 


78+1° 
80l?4 


20tl 
24+5 
311^ 



observed in this experiment are shown in Fig. [8] along 
with the predictions of the EFAX parameterization. The 
systematic errors associated with the target thickness, 
beam current measurements, and the widths of the mo- 
mentum slits have been included in the error bars and 
were 11% for the first data set, and approximately 2% 
for production data sets (8-13). For those isotopes that 
relied on transmission calculations, the weighted mean of 
all measured yields was used to obtain the "model-based" 
cross section (shown by solid diamonds in Fig. [8]) . The 
uncertainties in these cases included the statistical, the 
systematic and the transmission uncertainties. As can 
be seen in Fig. [51 the model-based cross sections are in 
good agreement with those produced by integrating the 
measured longitudinal momentum distributions as was 
shown above for the examples ^'^Cr and ^°K. 

The runs with the tungsten target employed strip- 
per foils (Be and C) placed directly behind the target 
which should be taken into account for the determina- 
tion of the cross sections. Such a correction can easily 
be made for the most neutron-rich isotopes (for example 
59-61 -J .^yi^ic]^ ^}jg ratio of the observed cross sec- 
tions R = CFw/o'Be exceeds «10. Cross sections with 
the tungsten target for isotopes with i? wl can be reli- 
ably obtained if the corresponding cross section with the 
beryllium target was measured with a small overall un- 
certainty (for example '^^'^^Sc). On the other hand, the 
division of the yields of ^''"^^Sc measured with the tung- 
sten target and a stripper between the two target materi- 
als was very uncertain. In addition, reliable lower limits 
could not be established for the two interesting isotopes 
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Qg (MeV) 



FIG. 9: (Color online) 
Cross sections for the 
production of neutron- 
rich nuclei with (a,c) 
odd atomic numbers 
and (b,d) even atomic 
numbers, with (a,b) a 
beryllium target 
and (c,d) tungsten tar- 
get, respectively. See 
text for explanation of 
Qg and the lines. The 
cross sections for ^^Ti at 
the center of the pro- 
posed new island of in- 
version [ll| are circled. 



^•^Al and ^^Mn due to large uncertainties in the transmis- 
sion. The cross sections for the production of some other 
isotopes (e.g., ^"'^Cl and ^^K) with the beryllium targets, 
needed to correct the measured cross sections with the 
tungsten target, had to be extrapolated. Finally, the 
cross sections needed for the correction of the data with 
the carbon stripper (9*'' data set in Table |TT| were calcu- 
lated by normalizing the beryllium cross sections by the 
ratio of the total reaction cross sections (Jc/crBe = 1-05 
from the EFAX parametrization. 



The compilation of results in Fig. [5] clearly shows a 
larger (sometimes exceeding a factor of ten) cross sec- 
tion with the tungsten target for the production of very 
neutron-rich isotopes of the elements 16 < ^ < 25 at 
this projectile energy. The geometrical factor (based on 
EFAX) is only 1.9. Similarly large ratios were recently ob- 
tained with a ^*Ca primary beam at a very similar beam 
energy Q . Models of nuclear reactions used for counting 
rate estimates, like the intranuclear-cascade plus evapo- 
ration model [2^ or abrasion-ablation in LISE++ [26] do 
not reproduce the low yields of exotic nuclei observed in 
this work. The predictions of the EFAX parameterization 
for reactions with beryllium, shown by the solid lines in 
Fig. [U reproduces the measured cross sections rather well 
in the region 13 < Z < 16. The data for the elements 
with Z > 16 generally fall below the predicted values but 
have a similar slope. 



C. Qg systematics 

The production cross sections for the most neutron- 
rich projectile fragments have been previously shown to 
have an exponential dependance on Qg (the difference in 




-70 -60 -50 -40 

Q (MeV) 



FIG. 10: (Color online) Production cross sections of neutron- 
rich calcium isotopes with a beryllium target as a function of 
Qg calculated with different mass models, see text. 
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mass-excess of the beam particle and the observed frag- 
ment) However, it is important to note that the 
actual masses of the very neutron-rich nuclei, needed for 
this calculation, have not been measured and are only 
available from models. There are many differences among 
mass models therefore deviations from the predicted yield 
using a given mass model might be able to identify miss- 
ing features in that mass model. For example, significant 
deviations of the real masses can be expected from the 
model predictions if a new island of inversion near ^^Ti 
(m (similar to the island of inversion observed near ■^^Na) 
is present but is omitted from the models. In a shell- 
model picture, the ground states of nuclei in the new 
island of inversion would be dominated by intruder con- 
figurations corresponding to neutron particle-hole excita- 
tions across the = 40 sub-shell gap into the 39/2 [IH- 
As a test of this suggestion, the cross sections for the 
production of all nuclei observed with the beryllium and 
tungsten targets arc shown in Fig. [5] where the abscissa, 
Qg, is the difference between the mass of the ground state 
of the projectile and the observed fragment. In this case 
the masses were taken from Ref. [Ifll- This empirical 
mass model (KTUY) provides an excellent systematiza- 
tion of the variation of the data, with the logarithm of 
the cross sections for each isotopic chain falling on an 
approximately straight line. 

The cross sections for each isotopic chain were fitted 
with the simple expression: 



<7(Z,A) = /(Z)cxp (Qg/T), 



(3) 



where T is an effective temperature extracted from the 
slope. Nearly all the data from the beryllium targets 
can be fitted with the single value of T = 1.8 MeV, but 
slightly better fits are obtained by finding the optimum 
T for each element. However, the heaviest isotopes of 
elements in the middle of the distribution {Z = 19, 20, 
21, and 22) appear to break away from the straight-line 
behavior. These heaviest four or five isotopes of these 
elements were found to have a shallower slope (shown by 
the solid lines) or enhanced cross sections relative to the 
model prediction. 

The systematic variation of the production cross sec- 
tions as a function of Qg was checked with several other 
well-know mass models and essentially the same be- 
havior was observed. Note that the differences among 
the mass models are small compared to the abscissa in 
Fig. [9] Specifically, models based on the Hartree-Fock- 
Bogohubov method HFB-8, HFB-9 ^ including the 
very recent HFB-17 version [2^, the finite range droplet 
model (FRDM) and the 2003 Atomic Mass Evalua- 
tion (AME2003) [3^ with shell crossing corrections that 
were developed in the LISE^"*" framework 31] were found 
to give similar results. The different behavior of the cal- 
cium isotopes with the beryllium target can be seen more 
clearly in Fig. [TOl The same effect is present in the same 
mass region with the tungsten target. 

The variation of the individual fitted values of the in- 
verse slope parameter, T, for products from both targets 




FIG. 11: (Color online) Values of the inverse slope parameter, 
T, from the best fit of Eq. |3]to the experimental cross sections 
in Fig. [9] shown as a function of atomic number. Open trian- 
gles connected by dotted lines are from the tungsten target, 
half-filled diamonds connected by dashed lines are from the 
beryllium target and the filled circles connected by solid lines 
are from the heaviest isotopes from the beryllium target in 
Fig. [9] 



is shown as a function of atomic number in Fig[TT] The 
inverse slopes of the cross sections from the beryllium 
target connected by dashed (solid) lines in Fig. [9] are in- 
dicated by the half-filled black diamonds and dashed lines 
(the filled circles and solid lines) in Fig. [11] Similar for 
the tungsten target the inverse slopes of the measured 
cross sections are shown by the open red triangles con- 
nected by dotted lines in Fig. [Til but the range of cross 
sections was more limited and only one slope parameter 
was fitted to their data. Fig. [S] clearly shows that there 
is a general increase in T for all of the heavy isotopes 
observed with Z = 19, 20, 21, and 22. 

One possible explanation for the shallower slope of the 
heaviest nuclei of these isotopes is that these nuclei are 
more bound (i.e., less negative Qg) than predicted by the 
mass models. A source of stronger binding can be defor- 
mation. In a shell-model framework, the wave functions 
of the ground and low-lying excited states of nuclei in the 
new island of inversion around ^■^Ti would be dominated 
by neutron particle-hole intruder excitations across the 
iV = 40 sub-shell gap, leading to deformation and shape 
coexistence. Thus, mass measurements of these nuclei, 
although difficult, would be especially interesting. 



V. SUMMARY 

The present study of the fragmentation of a ^^Ge beam 
at 132 MeV/u provided evidence for the production of 
fifteen previously unobserved neutron-rich isotopes p^ . 
The momentum distributions and cross sections for a 
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large number of neutron-rich nuclei produced by the Ge 
beam were measured by varying the target thickness 
in a two-stage fragment separator with a fixed momen- 
tum acceptance. The longitudinal momentum distribu- 
tions of 34 neutron-rich isotopes of the elements with 
13 < Z < 27 were compared to models that describe 
the shape and centroid of fragment momentum distribu- 
tions. New parameters for the semiempirical momentum 
distribution model [l9| based on the measured momenta 
were obtained. The production cross sections support 
the previous analysis on the basis of the Qg function but 
the most neutron-rich nuclei of elements with Z = 19 to 
22 are produced with an enhanced rate compared to the 
systematics. This trend was previously reported for frag- 
mentation on a beryllium target [l^ and was found to 



also be present with the tungsten target. The enhanced 
yields lie in the region that was previously predicted to be 
a new island of inversion in which an intruder shell-model 
state becomes the ground state. 
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APPENDIX A: IDENTIFICATION OF HEAVY 
IONS WITH ENERGY LOSSES 

The magnetic rigidity [Bp) of each particle was ob- 
tained by combining an angle measurement provided by 
two PPAC's in the dispersive plane of the S800 analy- 
sis beam line in combination with NMR measurements 
of the dipole fields: Bp — Bpo{l + S), where Bpo is the 
magnetic rigidity of the dipole at the central axis, and 
S = Ap/p is the fractional deviation of the particle from 
the central rigidity (for a constant charge q). The mass- 
to-charge ratio A/q is given by the expression: 



A/q 



Bp e 
uc 



(Al) 



where e is the elementary charge, u is the atomic mass 
unit, (3 is the velocity of the ion relative to the speed 
of light c, and 7 is the Lorentz factor. The velocity of 
the ions was determined from the time-of-flight (TOF). 
However, for the TOF measurement between the A1900 
focal-plane scintillator and the second PIN detector (see 
Table [T] and Fig. ^ , there were several materials in the 
path of the beam leading to different ion velocities along 
the beam line. In order to calculate the A/q ratio in 
this case with Ea lAU it was necessary to introduce a re- 
duced Bpa value that corresponds to the average velocity 
f^aver = L / (cT) , whcrc T is the TOF between two timing 
detectors. The time-of-flight value T can be written as 
T = Li/{c(3i), where N is the number of stages where 
the ion velocity changed due to the passage through ma- 
terial, and Li is the corresponding flight path at a given 
velocity. 

At intermediate beam energies when the ion's velocity 
is not significantly changed along the path (e.g. only 
thin materials are inserted in the beam) one can use the 
reduced magnetic rigidity as: 



N 



Bpa=kL/J2iL^/Bpo^ 



(A2) 



where k is the momentum deviation measured before the 
last fiight stage, and BpQi is the magnetic rigidity corre- 
sponding to the z'^'-stage on the optical axis. At these 
energies {f3 = 0.4) it is possible to omit the 7-factor 
in Eq. EI] because the ratio of {6 13 / (3) / {5^ / ^) = 5.2, 
whereas it is equal to 1.0 aX, (3 — 0.71. 

Calculations using Eq. IA2I for ^*Ca fragments (Data 
set 13, see Table |TT| with the TOF measured between 
the FP scintillator and the 2'"* PIN detector (see Ta- 
ble [l| were compared to LISE++ Monte Carlo simula- 
tions. A maximum difference of d(3 = 0.06% occurs for 
ions farthest from the central rigidity {Ap/p = 2.5%). 
This results in a mass resolution for ^^Ca of about 0.04 
units. 

The charge state (g) of the ion was evaluated with a re- 
lation based on the total kinetic energy (TKE), velocity, 
and magnetic rigidity: 



TKE 



u{j-l)iA/q) 



(A3) 



where TKE is calculated as a sum of the energy loss val- 
ues in each of the Si detectors. The velocity before the Si 
telescope was calculated from the reduced magnetic rigid- 
ity and the central magnetic rigidity of the last dipole in 
the beam line (Bpof): 



k Paver Bpof 



(A4) 



The atomic number was determined from the combi- 
nation of energy loss (AE) and velocity values according 
to the Bethe formula: 



AE 



1 , 5930 
-^log 



[MeV] " 1//3' 



(A5) 



